Bi-collinear antiferromagnetic order in the tetragonal a-FeTe 
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By the first-principles electronic structure calculations, we find that the ground state of PbO- 
type tetragonal a-FeTe is in a bi-collinear antiferromagnetic state, in which the Fe local moments 
(~ 2.5/is) are ordered ferromagnetically along a diagonal direction and antiferromagnetically along 
the other diagonal direction on the Fe square lattice. This bi-collinear order results from the interplay 
among the nearest, next nearest, and next next nearest neighbor superexchange interactions Ji, 
J2, and J3, mediated by Te 5p-band. In contrast, the ground state of a-FeSe is in the collinear 
antiferromagnetic order, similar as in LaFeAsO and BaFe2As2. 

PACS numbers: 74.25.Jb, 71.18.+y, 74.70.-b, 74.25.Ha, 71.20.-b 



The recent discovery of superconductivity in the lay- 
ered iron-based compounds attracts great interest, not 
only because the second highest superconductivity tem- 
peratures are achieved, but also because it is the first 
time that Fe atoms directly play an important role in 
superconductivity while the Fe atoms are convention- 
ally considered in disfavor of superconductivity. Un- 
til now there are four types of iron-based compounds 
that have been reported to show superconductivity when 
doping or under high pressures. They are classified as 
'1111' type (prototype: LaFeAsO [1]), '122' type (proto- 
type: BaFe 2 As 2 0), '111' type (prototype: LiFeAs H), 
and '11' type (prototype: tetragonal a-FeSe(Te) [4], [5[). 
The common feature shared by these compounds is that 
there are the robust tetrahedral layers where Fe atoms 
are tetragonally coordinated by pnictide or chalcogenide 
atoms and the superconduction pairing happens on these 
layers. 

It has been found that for both '1111' and '122' com- 
pounds there happens a tetragonal-orthorhombic struc- 
tural phase transition [6|. For undoped '1111', about 15K 
below this structural transition, a magnetic phase transi- 
tion takes place while for undoped '122' compounds, the 
structural and magnetic phase transitions occur at the 
same temperature. The magnetic order associated with 
the transition is a collinear antiferromagnetic order. The 
microscopic mechanisms underlying the structural tran- 
sition and antiferromagnetic transition and the relation- 
ship between these two transitions have been centered in 
hot debate. There are two very different microscopic pic- 
tures in understanding the phase transitions. The first 
one suggests that there are no local moments and the 
collinear antiferromagnetic order is entirely induced by 
the Fermi surface nesting which is also likely responsible 
for the structural transition due to breaking the four- 
fold rotational symmetry 0, @|. The second one suggests 
that As-bridged superexchange antiferromagnetic inter- 
actions between the nearest neighbor and next nearest 



neighbor Fe-Fe fluctuating local moments are the driving 
force upon the two transitions @, OH EH . 

In this letter, we report the electronic structure calcu- 
lations on teragonal a-FeTe and a-FeSe. We find that 
the ground state of a-FeTe is in a bi-collinear antifer- 
romagnetic order as schematically shown in Fig. 1(a), 
which differs from the collinear antiferromagnetic order. 
The underlying physics of this magnetic ordering can be 
described well in a model with the magnetic exchange in- 
teractions Ji, J2, and J3 respectively between the near- 
est, next nearest, and next next nearest neighbor medi- 
ated by Te-5p band (Fig. 1(a)). Hence, the fluctuating 
local moment picture works well in understanding the 
magnetism while the magnetic order can not be easily 
obtained from the Fermi Surface nesting picture. In the 
bi-collinear antiferromagnetic state, the bi-collinear anti- 
ferromagnetic ordering is along the diagonal direction of 
the Fe-Fe square lattice with double collinearity. If the 
Fe-Fe square lattice is divided into two square sublattices 
A and B shown in Fig. 1(a), the Fe moments on each 
Fe-Fe sublattice take its own collinear antiferromagnetic 
order. Our results are in excellent agreement with neu- 
tron experimental results [12L 1 13j|. In contrast, the ground 
state of tetragonal a-FeSe is a conventional collinear an- 
tiferromagnetic state, similar to LaFeAsO and BaFe2As2- 

In the calculations the plane wave basis method was 
used [14j. We employed the local (spin) density ap- 



proximation and the generalized gradient approxima- 
tion of Perdew-Burke-Ernzerhof [15[ for the exchange- 
correlation potentials. The ultrasoft pseudopotentials 
[l6| were used to model the electron-ion interactions. Af- 
ter the full convergence test, the kinetic energy cut-off 
and the charge density cut-off of the plane wave basis 
were chosen to be 800eV and 6400eV, respectively. The 
Gaussian broadening technique was used and a mesh of 
18x18x9 k-points were sampled for the Brillouin-zone in- 
tegration. In the calculations, the internal atomic coordi- 
nates within the cell were determined by the energy min- 
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FIG. 1: (Color online) Schematic top view of the FeTe (FeSe) 
layer. The ordered Fe spins are shown by red arrows. The 
small dashed square is an a x b (a = b) unit cell, (a) bi- 
collinear ordering with a 2a x b as the unit cell shown by solid 
rectangle; (b) collinear ordering with a \/~2a x y/2b as the unit 
cell shown by solid square. 



imization. Among all the four types of iron-based com- 
pounds, a-FeSe and a-FeTe possess the simplest struc- 
ture, made up only of two kinds of atoms Fe and Se or Te. 
They consist of continuous stacking of tetrahedral FeSe or 
FeTe layers along c-axis without any other kind of layers. 
a-FeSe and a-FeTe belong to a tetragonal family with 
PbO-type structure and space group of PA/nmm. The 
experimental tetragonal lattice parameters were adopted 
in our calculations 0, [H, [TH, 13, 17]. 

Our nonmagnetic calculations on both a-FeSe and 
a-FeTe exclude any possible structural distortions like 
Jahn- Teller effect. This means that a structural distor- 
tion on a-FeSe or a-FeTe if existing is very likely driven 
by magnetic forces. The calculated nonmagnetic elec- 
tronic energy band structures and Fermi surfaces are the 
same as those reported in Ref. jTsj. 

Now let's consider possible magnetic orders. If we di- 
vide the Fe-Fe square lattice into two square sublattices 
A and 5, there may be ferromagnetic, square antiferro- 
magnetic, or collinear antiferromagnetic order for the Fe 
moments on each of the two sublattices. The combination 
of these magnetic orders on the two sublattices yield the 
ferromagnetic, square antiferromagnetic, collinear anti- 
ferromagnetic, and bi-collinear antiferromagnetic orders 
on the Fe-Fe square lattice, as shown in Fig. 1. To 
calculate the bi-collinear antiferromagnetic order and 
collinear antiferomagnetic order, we adopt 2a x a x c and 
y/2a x \/2a x c unit cells, respectively shown in Fig. 1(a) 
and (b). For the other magnetic orders, the ax ax c unit 
cell is used. 

If the energy of the nonmagnetic state is set to zero, 
through the calculations we find that the energies of the 
ferromagnetic (Efe), square antiferromagnetic (Eafm), 
collinear antiferromagnetic (Ecol), and bi-collinear an- 
tiferromagnetic states (Ebi) are (0.183, -0.101, -0.0152, 
-0.089) eV/Fe for a-FeSe and (-0.0897, -0.0980, -0.156, 
-0.166) eV/Fe for a-FeTe, respectively. Thus the ground 
state of a-FeSe is a collinear-ordered antiferromagnetic 
state, similar to LaFeAsO and BaFe2As2- In contrast, 
the ground stat of a-FeTe is a bi-collinear ordered antifer- 



romagnetic state, being lower in energy by ~10 meV/Fe 
than its collinear ordered antiferromagnetic state. 

The magnetic moment located around each Fe atom 
is found to be about 2.2 ~ 2.6 /ig, similar to the cases 
of LaFeAsO and BaFe2As2, varying weakly in the above 
four magnetically ordered states. Since these local mo- 
ments are embedded in the environment of itinerant elec- 
trons, the moment of Fe ions is fluctuating. Nevertheless, 
we may here reasonably consider the spin of Fe ions is 
equal to 1. But the corresponding long range ordered 
moment should be smaller because the correlated effect, 
especially the strong competition between different anti- 
ferromagnetic states, has not been fully included in the 
density functional theory calculation. 

To quantify the magnetic interactions, we assume that 
the energy differences between these magnetic orderings 
are predominantly contributed from the interactions be- 
tween the Fe spins (S = 1), which can be modeled by the 
following frustrated Heisenberg model with the nearest, 
next-nearest, and next next nearest neighbor couplings 
Ji, J 2 , and J 3 , 

H = Ji^S z .S J + J 2 Si-Sj + Js ^'4 (!) 

(ij) «<ij»> 

whereas (ij), <C ij ^> and (((ij))) denote the summa- 
tion over the nearest, the next-nearest and the next-next- 
nearest neighbors, respectively. Since both a-FeSe and 
FeTe are semimetals, the above modeling may miss the 
certain contribution from the itinerant electrons. How- 
ever we believe that the above simply modeling captures 
the substantial physics on the magnetic structures. From 
the calculated energy data, we find that for a-FeTe J\ = 
2AmeV/S 2 , J 2 = 15.8meV/S 2 , and J 3 = 10AmeV/S 2 
while for a-FeSe J t = 71meV/S 2 , J 2 = 48meV/S 2 , and 
Js = 8.5meV/S 2 . Notice that the ferromagnetic state is 
quite low in energy in a-FeTe, which is reason why J\ is 
so small. 

It is well-known that for a J\-J 2 antiferromagnetic 
square lattice the frustration between J\ and J 2 will de- 
struct the Neel state and induce a collinear antiferro- 
magnetic order when J 2 > J\/2 as shown in Fig. 1(b), in 
which we notice that each site spin has two ferromagnetic 
connection neighbors and two antiferromagnetic connec- 
tion neighbors. Further inspection of Fig. 1 shows that 
J 2 and Js are the nearest and the next-nearest couplings 
for each of sublattice A and B whereas J\ is the connec- 
tion coupling between A and B. When Js > J 2 /2 and 
2J2 + 4 Js > Ji[l9[, a collinear antiferromagnetic order- 
ing will take place on each sublattice. It turns out that 
the bi-collinear antiferromagnetic ordering occurs on the 
square lattice as shown in Fig. 1(a), in which we no- 
tice again that each site spin has two ferromagnetic con- 
nection neighbors and two antiferromagnetic connection 
neighbors that doesn't cost more energy on J\ in compar- 
ison with the collinear antiferromagnetic order. Overall, 
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FIG. 2: (Color online) (a) Electronic band structure of the 
bi-collinear-ordered antiferromagnetic a-FeTe. The Fermi en- 
ergy is set to zero, (b) Brillouin zone, (c) Fermi surface. FX 
(FX) corresponds to the parallel(antiparallel)-aligned mo- 
ment line. 
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FIG. 3: (Color online) (a) Electronic band structure of the 
collinear-ordered antiferromagnetic a-FeSe. The Fermi en- 
ergy is set to zero, (b) Brillouin zone, (c) Fermi surface. FX 
(FX) corresponds to the parallel(antiparallel)-aligned mo- 
ment line. 



when Js > J2/2 and J 2 > Ji/2, by recursion we may thus 
expect that a bi-collinear antiferromagnetic state would 
be the ground state of the frustrated J1-J2-J3 Heisen- 
berg model (1), which is in good agreement with the de- 
rived Ji, J2, and J3 on a- FeTe and a-FeSe. Especially, 
we can show that J3 = J2/2 + (Ecol - Ebi)/^S 2 and 
J 2 = Ji/2 + (E AFM - E C ol)/±S 2 [2l|. The energy or- 
dering is E FM > Eafm > Ecol > E B i for a- FeTe and 
Efm > E B i > Eafm > E C ol for a- FeSe. By model 
(1) we may thus easily understand the complex magnetic 
structures displayed by a- FeTe and a- FeSe. 

When we further consider the possible spin-phonon in- 
teraction, it is expected that there would be further lat- 
tice relaxation, similar to spin-Peierls distortion, that is, 
being slightly longer along spin anti-parallel alignment to 
further favor antiferromagnetic energy and shorter along 
spin parallel alignment to further favor ferromagnetic en- 
ergy. Correspondingly, and 7 increase slightly respec- 
tively for the bi-collinear and collinear cases (Fig. 1(a) 
and (b)). As a matter of fact, we indeed find such extra 
small structural distortions for both a-FeTe and a-FeSe. 
For a-FeTe, increases to 92.03° with an extra energy 
gain of ~5 meV/Fe while for a-FeSe, 7 increases to 90.5° 
with an extra energy gain of ~2 meV/Fe. As results, the 
crystal unit cell of a-FeTe (a-FeSe) on FeTe (FeSe) layer 
deforms from a square to a rectangle (rhombus), as shown 
in Fig. 1. However, we find that this small lattice distor- 
tion affects weakly the electronic band structure and the 
Fe moments. 

Our calculations also show that for both a-FeTe and 
a-FeSe the Fe magnetic moments between the nearest 
neighbor layers FeTe (FeSe) prefer the anti-parallel align- 
ment but with a small energy gain less than 1 meV/Fe, 
similar to LaFeAsO. It is thus very likely that here the 
magnetic phase transition would happen below the struc- 
tural transition temperature. 

Fig. [2] shows the electronic band structure and the 
Fermi surface of a-FeTe in the bi-collinear antiferromag- 
netic order. There are three bands crossing the Fermi 



level which form three discrete parts of the Fermi sur- 
face. The Fermi surface contains a small hole-type pocket 
around R, two pieces of opened irregular hole-type sheets 
parallel to the plane T-Z-R-X, and two 'blood vessel'- 
like electron- type cylinders between F and X. From 
the volumes enclosed by these Fermi surface sheets, 
we find that the electron (hole) carrier density is 0.43 
electron/cell (0.41 hole/cell), namely, 2.38 x 10 21 /cra 3 
(2.26 x 10 21 /cm 3 ). In comparison with the nonmag- 
netic state, the carrier density is just reduced by half. 
The density of states at Ep is 1.98 state/ (eV Fe), which 
yields the electronic specific heat coefficients 7 = 4.665 
mJ/(K 2 * mol). 

The electronic band structure and Fermi surface of a- 
FeSe in collinear antiferromagnetic state is shown in Fig. 
[3l The Fermi surface consists of three sheets formed by 
two bands crossing the Fermi level. Specifically, one hole- 
type close sheet forms a hole pocket centered around Z 
and the other two electron-type cylinders are formed be- 
tween T and X, crossing the Fermi level. The volumes 
enclosed by these sheets give the electron (hole) carrier 
density 0.042 electrons/cell (0.028 holes/cell), namely 
2.68 x 10 20 /cra 3 (1.81 x 10 20 /cm 3 ). In comparison with 
the nonmagnetic state, the carrier density is reduced by 
one magnitude while in LaFeAsO the corresponding re- 
duction is more than hundred times. The density of 
states at Ep is 0.48 state/ (eV Fe), which yields the elec- 
tronic specific heat coefficients 7 = 1.129 mJ/(K 2 *mol). 

By projecting the density of states onto the five 3d 
orbit als of Fe, we find that the five up-spin orbit als are 
almost completely filled while the five down-spin orbitals 
are nearly uniformly filled by half. This indicates that 
the crystal field splitting imposed by Se or Te atoms is 
very weak and the Fe 3d-orbitals hybridize strongly with 
each other. As the Hund rule coupling is strong, this 
would lead to a large magnetic moment formed around 
each Fe atom, as found in our calculations. Moreover, 
Fig. 2] clearly shows that the band states constituted by 
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FIG. 4: (Color online) Top (left) and off-top (right) views 
of differential charge density for a-FeSe. The isosurface value 
was set to 0.03 e/ A 3 . Colors mapped on isosurfaces represent 
the relative height of data points in c-direction. 
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FIG. 5: (Color online) Total and orbital-resolved partial den- 
sity of states (spin-up part): a) bi-collinear antiferromagnetic 
a-FeTe; b) collinear antiferromagnetic a-FeSe. 

Fe 3d orbitals are very localized. Fig. [H plots the dif- 
ferential charge density distribution for a-FeSe (a-FeTe 
and LaFeAsO all have similar distribution [21]). We find 
that the most charge accumulations are surrounding Se 
(Te, As) atoms (larger pink balls), while the Fe atoms 
(smaller red balls) only have very small charge accumu- 
lation. Figured] (right) further shows a differential charge 
density pipeline connecting from a higher Se to a lower 
Se. For each Se atom, there are four charge accumula- 
tion pipelines connecting from it to its four adjacent Se 
atoms, suggesting an electron network mostly formed by 
delocalized electrons of Se through covalent bonding, in 
which the very localized Fe electron states are embedded. 

On the other hand, Fig. [5] shows that the band formed 
by Se 4p orbitals (also As 4p orbitals) is gapped at the 
Fermi energy while the band formed by Te 5p orbitals 
is partially filled. So there are itinerant 5p electrons at 
Fermi energy involved in mediating the exchange inter- 
actions in a-FeTe, likely, a long-range effect. This may 
explain why the coupling J3 is nearly zero for a-FeSe, 
LaFeAsO and BaFe2As2- 

It was reported 0, [H[ in experiment that at 105K the 
tetragonal a-FeSe compound experiences a similar struc- 
tural distortion as LaFeAsO, with a = b and 7 changing 
from 90° to 90.3° while at 45K the tetragonal a-FeTe 
compound also experiences a structural distortion but 
with a + b (a = 3.824 A, b = 3.854 A) and 7=90°. 



This is in excellent agreement with our calculations. Fur- 
thermore, by the neutron scattering Bao and coworkers 
(l2| found that an incommensurate antiferromagnetic or- 
der propagates along the diagonal direction of the Fe-Fe 
square lattice, but they also find this incommensurate or- 
dering is easily tunable with composition and locks into a 
commensurate order in the metallic phase. By our calcu- 
lation, it is clear that the magnetic order must be along 
the diagonal direction for a-FeTe. A physical picture sug- 
gested in our study is that the Fe moments are mediated 
by Te 5p-band and the origin J3 exchange coupling may 
be well induced through a RKKY-type mechanism. It is 
very likely that the excess interstitial Fe moments drive 
the bi-collinear order into the incommensurate order. No- 
tice that this diagonal antiferromagnetic order cannot be 
straightforwardly understood by the Fermi surface nest- 
ing picture. 

In conclusion, we have presented the electronic band 
structure of a-FeTe based on the first-principles density 
functional theory calculations. Our studies show that the 
compound a-FeTe is a quasi-2-dimensional bi-collinear 
antiferromagnetic semimetal with a magnetic moment of 
~ 2.5/iB around each Fe atoms. 
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